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XXVI. Further Note on the Propagation of Radio Waves over 
a Finitely Conducting Spherical Earth. 
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Natuurkundig Laboratorium der N. V. Philips’ Gloeilampenfabrieken, 
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Summary. 


In a former article (III.) it was pointed out that in the theory of the 
propagation of radio waves over a finitely conducting earth two different 
approximations for the many Bessel functions occurring in the problem 
can be introduced, viz., (a) the more simple tang. approximation (up till 
now used by most workers in this field), and (b) the less simple Hankel 
approximation, but which is required when the wave-length is small 
compared with the circumference of the earth. In ITI. the field on the 
surface of the earth due to an emitter on the ground was caleulated with 
an exclusive use of this better Hankel approximation, whereas in II., 
the superiority of the approximation not being fully realized yet at the 
time, curves for the field above the ground were calculated with the aid 
of both approximations. The results obtained there have now been 
fully recalculated with the exclusive use of the Hankel approximation 
resulting in the figs. 5 and 6, which therefore should replace the former 
figs. 19 and 20 of II. In these new figures the curves could be given 
ranging from the immediate proximity of the emitter up to well into the 
shadow region, due to the good fit existing in the transition region between 
the geometric optical formula (5) and the residue formula (46). Here 
the formerly introduced '' divergence factor" was numerically of great 
importance. 


$1. Introduction. 
In three former papers I., II., and III. +} we developed the theory of the 
propagation of radio waves round a spherical finitely conducting earth. 
Two different ways of attack were worked out :— 


(A) An extension of a series solution first given by G. N. Watson, 
which is best applicable in the ' shadow region " of the emitter 
and which, further, will be called the Residue-method, and 
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(B) an approximate geometric-optical ray expression obtained from 
the rigorous solution of the problem, which was best applicable 
in the lit region, and which will, further, be called the geometric- 
optical method. 


In III. we stressed the point that two different approximations to the 
Bessel functions occurring in the solution (A) could be made use of, viz., 
what we called (a) the tang. approximation (I. (16) and III. (10 a)), 
and (b) the Hankel approximation (I. (17), and III. (10 0)), of which, 
as pointed out in ITI., the latter is decidedly better for the practical 
condition (ka)? > 1, although the tang. approximation leads toa simpler 
treatment of the problem. We here use the same notations as in these 
former papers. 

In III. (figs. 3, 4, 5, and 6) we gave a complete set of practical curves 
for the field on the surface of the earth (h,—0) and when the emitter 
is also on the ground (h,—0). These data were all calculated with the 
aid of the Hankel approximation only and therefore constitute a very 
accurate result. 

As to the numerical results for an elevated emitter and (or) receiver 
(II., figs. 19 and 20), these were obtained partly still with the less accurate 
lang. approximation, the superiority of the Hankel approximation not 
being known yet at the time. We therefore have repeated these numerical 
investigations for the case of either or both the emitter and receiver at 
a height of h=100 metres and for o=10-%, e—4, but now with the 
exclusive use of the Hankel approximation. Moreover, we have taken 
the opportunity to apply our former geometric-optical solution (B) also 
numerically, and have obtained an excellent fit between the numerical 
results calculated according to the methods (A) and (B), so that we could 
considerably both extend and improve our former data II., figs. 19 and 20, 
which are therefore to be replaced by the present figures 5 and 6. 

In our former investigations we could show that, contrary to the general 
belief, the waves bend smoothly over the horizon and have still a con- 
siderable amplitude in the ‘ shadow region." The present extension 
of our former investigations fully confims this conclusion, and waves of 
the order of one centimetre are beyond the horizon even stronger than 
we thought before. The present investigation also fully confirms our 
former conclusion that a full cut-off at the horizon only occurs for waves 
of the order of one millimetre and shorter. 


§2. Recapitulation of Former Results. 


We showed that, as a consequence of the strong absorption in the earth, 
the total field I7,,, (even behind the horizon) could be approximated to 
by. the direct radiation I^ , , and the once reflected radiation O* , 9, where 
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it must be remembered that both these radiations penetrate by diffraction 
into the “ shadow region." This leads in case (A) to the series II. (77), 
viz. :— 
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On 
where (see II., (76)) the position of the orders n, is determined by 
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occurring in (1) may be called the height gain factors, which can be used 
to extend our formula III. (20 b) (which was meant for h, —h,—0) to 
ei 5X 
l 
(27a) 
which expression therefore constitutes our present basis as regards the 
method (A). 
On the other hand, the application of the geometric-optical treatment 
(B) (see IT. (91)) leads to 
gi aD gie, (Ra Ra) 
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being the sum of the incident ray and the reflected ray, the amplitude 
of the latter being given by the product of the “ divergence factor ” 
&ıı and the “ spherical reflexion coefficient " R,- 


$3. The Height-gain Factors fh) and f,(h;). 
In II. these factors f, were only explicitly given in the cases 9 —0 and 
8-—o where 
ib? 
7 (kya) 84/ kajk 3 —1 
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In the general case the fang. approximation leads to 
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and a similar expression for f (h). As explained above, in II. a combined 
use was made of the tang. and the, better, Hankel approximation, but, 
due to a resonance effect, instead of an improvement, a less accurate 
result was obtained, which caused in II. (99), (101), and (102) wrong 


factors to occur. Infact, Vy, and Vy, are to be replaced by xi 25, 


resp. Kd X2 — 27, and in the factors jv 57,2 the term 1 must be 


omitted*. Thus if we limit ourselves for the moment to the tang. approxi- 
mation only we obtain for h,~0 and h,=0 instead of (101) : 
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both being valid for y,2 and x.2>3 and where the former first factor 
2eP/Gk D) has been replaced by the slightly more accurate expression 
26/1316 3). 

But, as said before, we consider the Hankel approximation as the 
better one of the two, and with it we obtain from IT. (97) for f,(h,) : 
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As (see II., pp. 852 and 853) the H, of the nominator of (6 b) behaves 
differently when the argument is > or <1, ?. e., when A, OF «56225, 


* Mr. G. Millington, in a private communication, kindly drew our attention 
to this discrepancy. 
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we approximate these two cases separately. When h>56A? and for the 
case of ultra short waves, where | 5 | <1, (6 b) can further be simplified 
to 


; 3 2 
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where, for s > 4, the sum J',4--J' -13 may be approximated by 
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This formula (9 b) was derived from (6 b) by replacing H; in the nominator 
by the first term of its asymptotic development, whereas in the denomina- 
tor Hy), was first expressed in J,/, and Jı, and thereupon the first two 
terms of its Taylor development in à were retained. For h near 562?? 
more terms of the asymptotic development of H, are necessary. 

On the other hand, when 4-«563?5 we better go back to the original 
definition (3) of f,(h) which, making use of its differential equation (I. (19)) 
and of the above equation (2) for the zero's n,, together with the approxi- 
mation 
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being a simpler form for II. (98). Moreover, a further consideration 
shows that (10) leads to à minimum in the field as a function of height, 
which, for negligible displacement current, occurs at 


Cc 
Sees 1/2 P 
hmin — 7 Gema. s 


Qa 
confirming in a simple way a result of Eckersley and Millington *. 

The height-gain factors f,(4,) and f,(h,), appearing in (4b) and which 
belong to the Residue method (A), have thus been developed in a form 
suitable for numerical application. In the following paragraphs we shall 
consider in more detail the divergence factor «,, and the spherical reflexion 
coefficient R,,, both belonging to the geometric-optical approximation (B). 


* Phil. Trans. Roy. Soc. cexxxvii. p. 297 (1938). 


266 Drs. B. v. d. Pol and H. Bremmer on Radio 


$4. The Divergence Factor o4. 


The divergence factor «,, was developed in II. (91) from the general 
wave solution with the aid of the saddle-point method. It was explained 
physically as the additional divergence of a beam of rays due to the 
reflexion against a curved surface. The analytical expression found for 
this geometrical function was 


_ (RtR)Vkasinr — o(B,-- Ra) Vsin v, cos 7, 
Q4, == ———————————————— Áo SEE ES 
j tbra J sin 3. 2 .€08 7, COS Ty Vr sin $(Ear cos r,-- Rab cos 7) » 
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Fig. 1. 


Notations used in the geometric-optical formula (5). 


but as this factor is of great numerical importance in determining the 
field before the horizon we give here another, this time purely geometrical, 
interpretation of this function. 

Consider fig. 1, a narrow beam leaving the source Q and thereupon, 
after reflexion against the sphere at T, arriving at P. Let q be the actual 
normal cross-section of this beam at P; and let q' be the cross-section if the 


Waves over a Finitely Conducting Spherical Earth. 267 


reflexion had taken place against a plane surface at T, then we shall 


show that 
g 
CREE 
à q 


Obviously the cross-section a at P is given by 

—(R, -- R,)*.sin 7,.d7,d¢, 
where the angle 7, is as ea in fig. 1. The spatial angle of the beam 
on leaving Q is sin'r, dt, dd, where ¢ is measured at right angles to the 
plane of the paper. The cross-section q,, which the beam at P cuts out 
from a sphere S of radius r=OP and centre at O, is 


q, —1?.sin 3.d3d 4. 
Now, for fixed r, d .dn (see for definition of A II. p. 843) where 


n=k,b sin 7,. Further, we have dn=k,b cos 7,.d7,, from which it follows 
that 


and also q —4,.€08 T4. 


Having thus found q’ and q equation (12) at once gives æ, which is 
therefore seen to be identical with (11). Hence the divergence factor 
formerly derived with the aid of the saddle-point method can thus be 
fully explained geometrically. 

As to the numerical value of «,, this is unity on the lit part of the 
surface of the sphere, but it becomes zero on the line of sight and else- 
where varies between these two limits. 


$ 5. The Spherical Reflexion Coefficient Ry. 


As shown in (5), a numerical knowledge of the spherical reflexion co- 
efficient R,, is necessary in order to calculate the field in the lit part of 
space. The a for R,, found in H. a a) was 
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For a practical case, where a: position of the emitter and receiver 
are given, first the position of the reflexion point T (fig. 1) on the surface 
of the earth has to be determined, so that then, also, the angle 7,is known. 
Hence in (13) we must take n=k,a sin Ta, which can be written n=aà, 
where A=k, sin 7,. Our next task is to obtain approximations of the 


t= koe 
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t functions in (13). As further in the ¢‘ (ka) functions the argument 
ka is never near the order n=k,a sin 74, we have, as in I., p. 152, 
d M A— k, 
— lo gie x l ~ 2 
| + log (ett? (a)} NELLE 
whereas for the other ¢ functions, where the argument k,a may be near 


the order k, asin 7,, we must use the Hankel approximation according 
to I. (22). We thus obtain, instead of (13), 
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T'wo special cases are of interest :— 


(a) a> , leading to the plane problem. Here (14) gives at once the 
well-known Fresnel reflexion coefficient in the form as already 
considered in I. (21). 

(b) 7,—7/2 corresponding to grazing incidence. Here A=4,, and hence 
u=0 leading to the special case 

] —e-?/?1 AS 
peter AS’ 

605 17 (2/3) 

where A= -5 Hu) 220-9189. 

As for the ultra-short waves and average soil | 9| <1, we see 
from (15) that even for grazing incidence the spherical reflexion 
coefficient is very approximately equal to —1, which is also the 
value of the plane reflexion coefficient. 

The difference between the spherical and plane reflexion coefficient 
even at grazing incidence thus being very small, we may, for the con- 
ditions mentioned, everywhere in (5) replace the spherical by the ordinary 
plane reflexion coefficient. 'This, however, would not be allowed if 
longer waves or a reflexion against soil of a better conductivity were 
investigated. Æ. g., if the conductivity were infinite, (9—-oo ), (15) would 
give for the value of the spherical reflexion coefficient at grazing incidence 
€-*"? instead of +1 as for the plane case. 


(15) 


§ 6. Numerical Results. 
Above we have given all the details for a numerical evaluation of cases 
of practical importance, all deductions now being exclusively based on 
the more accurate Hankel approximation. 
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With the aid thereof we have again with great care constructed two 
sets of curves giving the attenuation (as compared with the free radiation) 
of the waves as a function of distance, viz. : 


(a) for h,=100 metres, h,=100 metres, 
(b) for h,=100 metres, ka —0, 


the conductivity c and dielectric constant e being taken o —101-? e.m.u., 
«=-4 representing average soil. The range of distance covered is from 
0 to well into the “ shadow region "; the curves have been calculated 
for the wave-lengths A=7 m., 70 cm., 7 cem., and 7 mm. 

As already stated above, the parts of the curves for the lit region were 
calculated with the geometric-optical formula (5), the divergence factor 
o4; being duly taken into account. On the other hand, the residue formula 
(4 b), which is principally applicable anywhere, is most suited beyond and 
near the horizon, and a range of overlap with (5) exists, showing a very 
good numerical agreement. Thus (45) and (5) together enable us to 
draw uninterrupted field curves for distances ranging from close to the 
emitter up to far into the shadow region. 

Before, however, describing the final curves we give three additional 
figures 2, 3, and 4, which show some detailed results of the calculation. 
Consider first fig. 2, which has been drawn for o=10-*8, e=4, A=0-7 metre, 
h,=h,=100 metres. In this figure, therefore, the emitter is at a height 
of 100 metres, whereas the receiver is thought to be all the time at the 
same height, the linear distance D being varied. Up to the shadow limit 
maxima and minima occur due to the interference of the direct and 
reflected radiation as represented by the first and second term of (5). 
Close to the emitter the ratio [/IT,, is obviously unity, but soon we reach 
the distance where it can be said that the interference maxima lie on 
the curve C, given by | JZ/II,, | —1--«,,| Ry, | and the minima on the 
curve C,, viz. | H/I | =1—%1 | Ry |. It is seen from fig. 2 that the 
distance 2a,, | R4, | between C, and C, shows a minimum at D —0-4 km., 
this being due to the Brewster minimum in the reflexion. It may be 
remarked that beyond this minimum the effect of the divergence factor «,, 
becomes more and more of importance ; if we had wrongly considered the 
reflexion to take place against a flat earth, the interference minimum at 
D —13 km., e. g., would have become 0-066 instead of 0-096. 

When we approach the shadow limit and penetrate into it the residue 
series formula (45) must be used. Therefore in fig. 2 all the points 
caleulated from the geometric-optical formula (5) are shown as black 
circles, whereas those obtained from the residue series as open circles. 
The numbers placed near these open circles represent the number of 
terms of the residue series which were necessary in order to obtain an 
accuracy of one per cent. The figure shows very clearly how this number 
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Two examples showing the evaluation of the residue series as a vectorial 
addition of its successive terms. The top drawing refers to a point between 
the last interference maximum and minimum, whereas the lower figure 
refers to a point between the last maximum and the shadow limit. 
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of terms rapidly increases from one well beyond the shadow limit, up to 
10 to 17 near the first interference maximum. Moreover, it demonstrates 
the excellent fit obtained in the overlapping range which occurs near 
this first interference maximum. 

Turning next to fig. 3 we there represent for two special cases the 
vectorial addition of the several complex terms of the residue series, 
viz., for o —10-13, e=4, A—0-7 m., h,—h4—100 m., and the distances 
D=29-75 and 46.28 km. The numbers placed near the corners of the 
broken spirals indicate the summation of 1, 2, .... terms of the series. 


Fig. 4. 


An example showing what becomes of the Cornu spiral of optics for 
the diffraction of radio waves round an absorbing sphere. 


In order to obtain with sufficient accuracy the final vector, indicated by 
the dotted line from the origin, thirteen terms were necessary in the top 
curve, whereas only four terms were sufficient in the case of the lower 
figure, which refers to a bigger distance. 

Coming next to the rather complicated fig. 4, which was calculated 
for o —10-33, e=4, A=7 m., h,=h,=100 m., this drawing represents what 
becomes of the well-known Cornu spiral in our special diffraction problem 
with absorption. The numbers 3 km., 5 km., 10 km..... etc. placed 
near the curve indicate the distance of the receiver from the emitter, 
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and the vector drawn from the origin O to any point on the spiral represents, 
both in amplitude and phase, the ratio of the actual field at the corre- 
sponding distance to the field which would obtain there in the absence of 
the diffracting earth. If we follow the spiral from the centre (1-072), 
which corresponds to the receiver near the emitter, it first spirals out and. 
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The attenuation factor for short waves when both emitter and receiver 
are at a height of 100 m. and for propagation over average soil. 


then shrinks again to a minimum value at (B) corresponding to the 
Brewster minimum, as fully explainedwith fig. 2; next it expands again 
(the radius corresponding all the time to «,, | Ri | ) till it bends over 
at the point S, where the shadow limit is reached, finally to approach 
the origin O again in a spiral way, but with inversed rotation, the details 


of this latter part being given in the left top drawing. 
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Having thus explained some details of our general numerical pro- 
cedure we give in figs. 5 and 6 the final results. Both figures have again 
been calculated for c=10-1, e=4, A=7 m., 0-7 m., 7 cm., and 7 mm. 
Fig. 5 relates to the case where both emitter and receiver are at 100 metres 
height (h,=h,=100 m.), whereas in fig. 6 the emitter is supposed at 
h,=100 metres, the receiver being on the ground (5,—0), where it may 
be remarked that, due to the principle of reciprocity, the last figure also 
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The attenuation factor for short waves of the field at the surface of the earth 
when the emitter is at a height of 100 m. and for propagation over average 
soil. 


refers to the case where the emitter and receiver are interchanged. Thus 
the present more accurate fig. 5 has to replace fig. 20 of IT. and the present 
fig. 6, fig. 19 of II. Fig. 5 clearly shows again that a wave A=7 metres 
penetrates well into the optical shadow region, no discontinuity whatever 
occurring at the shadow limit. As found before, only waves of a length 
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of the order of one centimetre or shorter exhibit a rapid cut-off when 
we pass into this shadow region. 

The general behaviour of the curves of fig. 5 can be described as 
follows :—-The final decrease for short waves starts already at the last 
interference maxima (the one nearest the shadow limit), the slope of the 
curve varying further only very little when we enter into the shadow region. 
It should be observed that this last maximum occurs the nearer the shadow 
limit the shorter the wave-length, it lying on this limit when A—0, for 
which case an infinite sharp cut-off occurs. 

Turning now to fig. 6, this has been drawn for the same constants 
of the earth and wave-lengths, only h,=100 m. and A,=0. Here 
obviously no interference maxima and minima occur, the curves being 
of a fully monotonic character with a complete absence again of any 
discontinuity at the horizon. Before the horizon the geometric-optical 
formula (5) was again used up to the transition region, after which the 
residue series (4 b) was applied. Moreover, in this case, where the receiver 
is on the ground, we have for all distances two simplifications in the 
general formula (5), viz., «,,=1, and D=R; +R». 

It is further seen that before the horizon the curves asymptotically 
approach the limiting curve for A—0, which latter case could be very 
simply treated with the aid of (5), because then the conduction current 
is negligible compared with the dielectric displacement current (k,?—>ew?/c?), 
leading quite simply to 

II 2e COS T3 
Hy > € COS Tot Ve—sin?r, ' 
where 7, is the angle of incidence of the direct ray at the receiving point. 

In conclusion, the feeling may be expressed that, as far as we are aware, 
the new curves given here, and which are again based on the hypothesis 
of à homogeneous atmosphere, constitute the most accurate and detailed 
theoretical results so far published because they were derived with one 
single approximation to the Bessel functions, viz., what we have called the 
Hankel approximation. 


Eindhoven. 
3lst October, 1938. 


